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We point out that the minimal seesaw model can provide a natural framework to accommodate 
tiny neutrino masses, while its experimental testability and notable predictiveness are still main- 
tained. This possibility is based on the observation that two heavy right-handed Majorana neutrinos 
in the minimal seesaw model may naturally emerge as a pseudo-Dirac fermion. In a specific scenario, 
we show that the tri-bimaximal neutrino mixing can be produced, and only the inverted neutrino 
mass hierarchy is allowed. The low-energy phenomena, including non-unitarity effects in neutrino 
oscillations, neutrinoless double-beta decays and rare lepton-flavor-violating decays of charged lep- 
tons la Ip^, have been explored. The collider signatures of the heavy singlet neutrino are also 
briefly discussed. 

PACS numbers: 11.30.Fs, 12.15.Ff, 14.60.St, 14.60.Pq 



I ■ 



On 
O 



X 



I. INTRODUCTION 



Recently, a lot of attention has been focused on the 
experimental testability of seesaw models for neutrino 
masses at the CERN Large Hadron Collider (LHC). In 
the typical extension of the standard model (SM), three 
right-handed neutrinos are introduced and assigned large 
Majorana masses P, H, [1, Q. In this case, the SM 
neutrinos acquire tiny Majorana masses via the type- 
I seesaw mechanism, i.e. the effective mass matrix of 
light neutrinos is given by w — M^Mp ^M^. Here 
the Dirac mass is naturally around the electroweak scale 



~ ^EW = -'^00 GeV and the heavy Majorana neu- 
trino masses are extremely large ~ ^O^" '^'' GeV. 
However, this type-I seesaw model suffers from the lack of 
testability, because right-handed neutrinos are too heavy 
to be produced in current collider experiments. 

In order for the typc-I seesaw model to be testable, 
we have to implement the structural cancellation condi- 
tion MpMp ^Afp « 0, which can lead to sub-eV neu- 
trino masses but keep heavy Majorana neutrino masses 
as low as several hundred GeV [H, Q. However, the fine- 
tunings of the structures of Mj^ and are unavoid- 
able, which seems unnatural. It has been pointed out 
that the heavy pseudo-Dirac neutrinos may account for 
tiny neutrino masses, such as in the so-called inverse see- 
saw models, which can be tested in collider experiments, 
and are free of the fine-tuning problem 0]. Neverthe- 
less, more singlet fermions should be added in order to 
form pseudo-Dirac fermions together with singlet right- 
handed neutrinos. For instance, in the minimal version of 
inverse seesaw model, we must have four singlet fermions 



to guarantee two massive light neutrinos [sj. Even more 
singlet fermions are required in the realization of multiple 
seesaw mechanisms which are the direct generations 
of the type-I and the inverse seesaw models. 

We propose that the minimal type-I seesaw model 
(MSM) with only two right-handed neutrinos may be 
the most natural candidate for realistic and testable neu- 
trino mass models at the TcV scale. The reason is simply 
that these two right-handed neutrinos themselves can be 
combined together to make a Dirac fermion in the U{1) 
symmetry limit. They may also be embedded into a two- 
dimensional representation of the discrete flavor symme- 
try groups, such as the permutation group S4 [HI. 
The soft symmetry-breaking terms then give rise to 
tiny neutrino masses, while the heavy pseudo-Dirac neu- 
trino provides us with rich low-energy phenomena, e.g. 
the non-unitarity effects in neutrino oscillations and the 
lepton-flavor-violating decays of charged leptons. Fur- 
thermore, the tri-lepton signals pp — *■ £^i^i^i'{D) -I- jets 
of the pseudo-Dirac neutrino can be discovered at the 
LHC. Due to the minimal number of model parameters, 
the observables at low energies and in the collider experi- 
ments are intimately correlated with each other, and then 
serve as a cross test of our minimal TeV seesaw model. 

The remaining part of this work is organized as follows. 
In Sccini we first present the structure of our model. The 
neutrino mass spectra and neutrino mixing patterns are 
discussed in detail in Sec. IIIII The implications for low- 
energy phenomena and possible collider signatures are 
addressed in Sec. IIVI Finally, a brief summary is given 
in Sec. El 
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In the minimal seesaw model [lH, , we extend the 
SM by introducing two heavy right-handed neutrinos. 
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which are singlets under the SM gauge group. The La- 
grangian relevant for neutrino masses is 



where = [v^ 



in^M^nR + h.c. , (1) 



•l)^ and 



''iR' "'2r)"^ stand 



for left- and right-handed neutrinos, respectively. To ef- 
fectively suppress the light neutrino masses while keeping 
heavy ones around the TeV scale, we additionally impose 
a global U{1) symmetry on the generic Lagrangian, under 
which the charges of the SM lepton doublets are opposite 
to n2j^ but equal to Tixp . With the help of such a Icpton- 
number-like symmetry, the mass matrices il/^ and il/^ 
in Eq. ([T]) take very simple forms 



Ve 





M 
M 



(2) 



with V ~ 174 GeV being the vacuum expectation value 
of the Higgs field. In the flavor basis (1/^, n^^), the overall 
5x5 neutrino mass matrix reads 



M = [ \ ' 



(3) 



Note that the rank of A^^ is two, so three light neu- 
trinos are massless and two heavy Majorana neutrinos 
are degenerate in mass, as a consequence of the addi- 
tional global symmetry. This can be easily verified by 
noting that the mass matrices in Eq. ^ satisfy the rela- 
tion Afj-iMj^^Mp = 0, which implies that light neutrino 
masses are vanishing to all orders. Therefore, the masses 
of heavy right-handed neutrinos have nothing to do with 
light neutrinos and can be located at a relatively low 
scale, e.g., around the TeV scale. On the other hand, the 
symmetric matrix Mp can be diagonalized via an orthog- 
onal transformation V^M^V^ = Diag{— A/, A/} with 



1 
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1 1 

-1 1 



(4) 



Thus the mass eigenstates of heavy right-handed neu- 
trinos Pi and P2 possess identical masses but opposite 
CP parities, and they constitute a four-component Dirac 
particle P ^ {P^ + P^I\/2 with mass being M |l4l | . 

In order to accommodate light neutrino masses, one 
can add small perturbations \,o M.^, which softly break 
the global V(X) symmetry. There are in principle four 
classes of soft-breaking perturbations to A^^, and the 
most general form of the perturbed neutrino mass matrix 
is given by 



1 ^ee 


















vye 




vyr 







vye 



vyr 

M 



M 



(5) 



All the above perturbation terms k^^, e^, fj,' and fi (for 
a, (3 = e, fi, t) break the lepton number conservation, and 



hence bring in neutrino masses proportional to the cor- 
responding couplings. Some comments on the possible 
origins of these perturbations are in order: 

(i) The K term corresponds to a purely Majorana mass 
term of light neutrinos, which can be realized in a more 
complicated theory with additional contributions to neu- 
trino masses. A typical example is the type-(I-l-II) see- 
saw model [H, [la, where an SU{2) triplet Higgs 
with mass much larger than the electroweak scale is in- 
volved. At lower-energy scales, the decoupling of the 
triplet Higgs will result in light neutrino masses together 
with non-standard interactions through the tree-level ex- 
change of the neutral scalar [l^ . Another possibility is to 
incorporate extra SM singlet or triplet fermions, which 
give birth to a neutrino mass operator similar to that 
in the type-I or type-Ill seesaw models [l^. Although 
feasible, the above mechanisms are always pestered with 
too many parameters, which render the models neither 
predictive nor economical. 

(ii) The /x term in Eq. ([S]) is a bare Majorana mass 
insertion violating the lepton number by two units, which 
is also realized in the inverse seesaw framework 0, ■ 
In the presence of the fi term, the light neutrino mass 
matrix can be obtained from the inverse seesaw formula 



M„ 



(6) 



Hence the smallness of neutrino masses is attributed to 
both the small ^ parameter and the ratio M-^/M. How- 
ever, as pointed out in Ref. [1], at least two pairs of sin- 
glet heavy neutrinos are required in order to enhance the 
rank of A^j^ from two to four. One can also see this point 
from Eq. ([6]) that the rank of Af^ is exactly one, which 
definitely comes into conflict with the observed two mass- 
squared differences in neutrino oscillation experiments. 

(iii) The fi' term in Eq. ^ does not contribute to neu- 
trino masses at the tree level. However, it may radiatively 
generate neutrino masses via one-loop diagrams involving 
right-handed neutrinos and gauge bosons [5| . In addition, 
due to the corrections induced by the n' term, the masses 
of Pi and P2 are not exactly equal, and thus the small 
mass splitting between heavy neutrinos could natu rally 



make the resonant leptogenesis mechanism feasible [21 
In analogy with the /x-term corrections, the drawback is 
that only one light neutrino may acquire mass, and hence 
the fj,' term is not phenomenologically adequate. 

(iv) The £ term softly violates the extra U{1) symme- 
try but enhances the rank of A^^, from two to four, which 
is required by the neutrino oscillations. Such a coupling 
could be easily realized in grand unified theories, e.g. in 
the supersymmetric 50(10) model with a very low B — L 
scale [23l. As we will show later in this case, neutrino 
masses are naturally tiny, since they are proportional to 
£, and are further suppressed by the mass ratio A/j-,/A/. 
In the following, we will only concentrate on this par- 
ticularly interesting pattern of neutrino mass generation, 
and figure out the phcnomcnological consequences in de- 
tail. The overall neutrino matrix can be obtained from 
Eq. (O by setting all k, fi and ^' to zero. 
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Without loss of generality, one can always redefine the 
lepton fields so as to remove the corresponding phases 
of and M , leaving only complex. Therefore, we 
shall assume and M to be real throughout the follow- 
ing discussions. The total neutrino matrix M^^ can be 
diagonalized by the unitary transformation 



It is obvious that the neutrino mixing matrix N ap- 
pearing in the charged current is non-unitary. As we 
shall show, the mixing between light and heavy neutrinos 
will bring in several significant phenomenological conse- 
quences, in particular, when the scale of the heavy neu- 
trino masses is accessible to future colliders. 



M,V* ^ M, = Bmg{m^,m^,m^,-M,M} , (7) 

where is a 5 x 5 unitary matrix, and (for i = 1,2,3) 
are masses of three light neutrinos. In the leading-order 
approximaiton, the effective mass matrix of light neutri- 
nos is given by the type-I seesaw formula 



M„ 



(8) 



where e = (e^, e^, e^)'^ and F = u){y^^ y^, y^.)'^ with ui = 
v/M. In general, M^^ can be diagonalized by a 3 x 3 
unitary matrix as U'^M^U* — = Diagjm]^, mj, rrig}, 
and U is usually parametrized in the standard form 



^P5Rl3i0l3)Ps''Rl2{0l2)PM , (9) 



where Rij correspond to the elementary rotations in the 
ij = 23, 13, and 12 planes, are the rotation angles, 
while Pg = Diag{l, l,e'*} and = Diagjl.e'P, 1} con- 
tain the Dirac and Majorana CP-violating phases, re- 
spectively. Note that only one Majorana phase is needed 
to parametrize U since one light neutrino is massless, 
which is the salient feature of MSM. The phases in 



P4, ^ Diag(e"^ 



are usually rotated away in 



the SM context but must be kept in the current model. 
It should be noticed that U is not exactly the matrix 
governing neutrino oscillations, even if we choose a ba- 
sis where the charged-lepton mass matrix is diagonal. To 
clarify this point, we turn back to the 5x5 unitary matrix 
V in Eq. which can be rewritten in a block form 



V 



N3X3 
S2x3 



R3X2 
^2x2 



(10) 



The approximate expression of each block can be found 
in Ref. [2^, and the relevant ones are 

Nc^(^l-^FF^^U, Rc,QV^, (11) 

with Q = (0, F) being a 3 x 2 matrix. The flavor eigen- 
states of neutrinos are then the superpositions of light 
neutrino mass eigenstates v^^^ = (i^iLi f^2Li ^3l) ^^^'^ 
heavy one P. More specifically, the flavor eigenstates of 
light neutrinos can be expressed as — Nv^^ + FP'^, 
which indicates that P"^ mixes with left-handed neutri- 
nos and enters into the weak interactions after the elec- 
troweak symmetry breaking. Therefore, the leptonic 
charged-current interactions in the mass basis read 



FP") W- + h.c. . (12) 



III. NEUTRINO MASSES AND MIXING 

First, we should consider the neutrino mass spectra 
and flavor mixing matrix. It is straightforward to obtain 
the light neutrino mass matrix from Eq. ^ 



M^=uj\ e^y^ + e^y^ 2e^y^ e^y^ + e^y^ 
yE^y^+e^y^ e^y^+e^y^ 2e^y^ 



,(13) 



where the irrelevant minus sign has been omitted. Since 
the experimental data on neutrino oscillations suggest 
the tri-bimaximal mixing pattern [2^ . [2^ , a /x-r symme- 
try is particularly favorable in constructing the neutrino 
mass matrix. To this end, we assume here that the rela- 
tions £^ = Et and y^i = j/t hold. Defining A = 2uje^y^, 



B = ^{SeVt, + e^ye) and C 
Eq. Uni) as follows 



M., = 



2(jj£^y^, one can rewrite 




(14) 



which can be further put into a 2 x 2 block form by a 
maximal rotation 

/ A V2B 0\ 
< = ^23(7)^^^3(7) ^\V2B 2C . (15) 
^ * V 0/ 

We have found that only the inverted mass hierarchy 
> ™i > ™3 = is compatible with a maxi- 
mal mixing pattern in the 2^3 sector. The left-up 
block matrix in Eq. (|15p . denoted as M^, is a general 
2x2 symmetric matrix, which can be diagonalized as 



I/qM^Uq = Diag{ 771^,7712} with 



e'^^i 
e''^'^ 



(16) 



Here Cg = cos 9 and Sg = sin 6* have been defined. After a 
lengthy but straightforward calculation, one can obtain 



X /2(X'^ + Y"^) 

tan <^ = - , tan 26 = ^ _ 4|g|2 ' (1^) 



where 



X = Im(S)Rc(.4 - 2C) - Re(B)Im(A - 2C) , 

Y = lm{B)lm{A + 2C) +Rc{B)Rc{A + 2C) , (18) 
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together with the mixing angles ^13 = and 6*23 = 45° de- 
fined in Eq. ^ . The non- vanishing neutrino mass eigen- 
values are given by 



mi 

1712 



Asl + V2Bsie-'"^ + 2Ccle-^'"^ 



(19) 



In addition, the phase difference A?/' = ?/'2 — "01 reads 



A-0 = 2 arg 



Asl + V2Bsle-'"t> + 2Cc^e-2i</' 



Ad 



(20) 




In comparison with the standard parametrization in 
Eq. ([9]), we can observe that p = Aip — is just the phys- 
ical Majorana phase. Once the identity VA^ + Y^ = 
2(|^P — 4|Cp) is fulfilled, the tri-bimaximal mixing pat- 
tern with 6*12 = 6* ~ 35.3°, 6I23 = 45° and 9^^ = is 
reproduced in the leading order. The vanishing ^13 can 
be regarded as a consequence of the exact ^ r sym- 
metry. If we relax this assumption, i.e. 7^ and 
(or) 7^ j/^, both non- vanishing mixing angle ^13 and 
CP-violating phase S can be accommodated. 

As we discussed above, one peculiar feature in our 
model is that the mixing matrix of light neutrinos N is 
no longer unitary. Adopting the parametrization of non- 
unitary leptonic mixing matrix = (1 — 7/)[/ in Ref. [26| . 
we can see from Eq. pT|) that 



(21) 



The /i <-> T symmetric feature of the effective neutrino 
mass matrix M^^ indicates 77^^ ^ rj^^. Since the most 
stringent experimental constraints come from the lepton- 
fiavor-violating (LFV) decay /i ej, severe bounds on 
T^g^ and 1]^^ are expected. 

It is worthwhile to note that the ratios e^/e^ and y^/y^ 
are of great importance in our model. For a given mass 
M (or equivalcntly oj = v/M) and also the mass scale of 
light neutrinos wi/g|eg|, all the other observables are de- 
termined by these ratios. One can always appropriately 
rescale and y^ to get the desired masses m^, keep- 
ing the ratios s^je^ and y^/yf, unchanged. Therefore, 
we shall mainly concentrate on the relative ratios among 
model parameters in the following discussions. 

In FIG. [1] we show the allowed parameter space of the 
model within 3(t C.L. In our numerical analysis, we take 
M = \ TeV for example, as well as the values of neutrino 
masses and mixing angles from Ref. [13] ■ We also include 
the experimental constraints on the non-unitary effects 
coming from universality tests of weak interactions, rare 
leptonic decays, invisible width of the Z-boson and neu- 
trino oscillation data [2^. Technologically, we randomly 
choose the values of (e^, Z/ei 2/^) ^'^'^ their correspond- 
ing phases, while the data sample reproducing all the 
neutrino masses and mixing angles within 3<t confidence 




2| / I 4 



180 



135 




FIG. 1: The parameter space of the neutrino mixing angle 
O12, the Majorana phase p, and the ratios (e^/Ce, T^^r/^yep)- 



ranges will be kept. The absolute scales of and y^ 
(for a = e,/i) are also checked to be consistent with the 
experimental bounds on 77. Some comments are in order: 

• From the uppermost plot, one can see that the ratio 
IVfir/Vefil = IVfi/yJ is Strictly constrained by 6*12 at 
large values, and the maximal value of the ratio is 
close to 5. In this case, the current bound \r]e^\ < 
6 X 10~^ suggests a rather stringent bound \i]f^r \ ^ 



5 



3 X 10~^. However, as pointed out in Ref. [1^, in 
the case with M < Aew but above a few GeV, 
the severe constraints from Z and W decays do not 
apply since the unitarity is restored. In this special 
case the upper bound |77e^| < 9 x 10~^ indicates 
that there is no general constraint on It^^^tI- 

The correlation between |£;^/£e| and \i]fiT/Vefj.\ is 
illustrated in the middle plot. The approximate re- 



lation |e^/£e| X \riefj,/r]fj 



0.5 reflects the fact that 



\A\ and 2\C\ in Eq. are comparable in magni- 
tude to ensure the experimental result rrii ~ mj. 
To generate a large but non-maximal mixing angle 
9^2, the relation \B\ <^ \A\ must hold, which is only 
possible if the relative phase between A and C (or 
equivalently and e^) is around 180°. However, 
the exact identity \A\ = 2\C\ is forbidden due to the 
mixing angle 9 in Eq. p7|) , which exhibits the strict 
experimental constraints on model parameters. 

• In the last plot, we show the allowed ranges of the 
Majorana phase p versus the ratio {ij^r/Venl- As we 
have discussed, the requirement A ~ — 2C also sets 
a strong constraint on the relative phase between 
Ee and e^. For a sizable ry^r, which is phcnomcno- 
logically interesting, p ~ 90° can be read from the 
plot. For intermediate values of 77^^, wc can see 
that p deviates significantly from 90°. In principle, 
the whole range [0, 180°] of p is allowed. 

As a consequence of the minimal number of model pa- 
rameters, the mixing angle 9^2, the Majorana phase p 
and the non- unitarity parameters \r]^p\ are now related 
to one another. Furthermore, the introduction of non- 
vanishing 0]^3 and CP-violating phase is straightforward 
by relaxing the p-T symmetry. 



IV. FURTHER DISCUSSIONS 

Now we shall address the possibly intriguing low- 
energy phenomena of our model, as well as the signals 
of the heavy neutrinos at the LHC. 

a. Rare LFV decays mediated by the heavy neutrino. 
The heavy neutrinos entering into the charged-current 
interaction will contribute to the LFV decays of charged 
leptons. For instance, the decay channel £^ ipj can 
be mediated by heavy neutrinos [2^. In the standard 
type-I seesaw scenario or MSM, one has approximately 
(Mj3/M)2 = O(mM-i), and therefore BR(^„ ^ ipj) cx 
0{m^) indicates a strong suppression of LFV decay rates, 
where m and M denote respectively the scales of light 
and heavy neutrino masses. However, in our model, one 
can have sizeable Mj^/M without facing the difhculty of 
neutrino mass generation since they are suppressed by 
the perturbations e^. Thus appreciable LFV rates could 
be obtained even for strictly masslcss light neutrinos [sO] ■ 

h. Neutrinoless double beta decays (Qv2[3). The 
light Majorana neutrinos contribute to Qv2[3 decays, 
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FIG. 2: The allowed region of the effective neutrino mass 
(m)^^ versus the ratio |r?^^/r/e^|. 



which serve as the unique tool to discriminate between 
Majorana and Dirac nature of massive neutrinos. The 
relevant quantity is the effective neutrino mass {m) = 

l^^i^ei + "^2^e2 + "^s^esli which in our model with 
7713 = can be evaluated as {rn) = \miU^i + rn2U^2\- 



Because of > "^1 VT^"^3il ^ 0-05 eV, the two 
terms in the {m)^^ are comparable in magnitude. In 
this case, the Majorana phase plays a key role in de- 
termining {m)pp. For example, if p is far away from 
90°, the contributions from these two terms should be 
added constructively, and one can then expect a large 



value 



I HP 



0.05 eV. In FIG. [2] we have shown the 



allowed region of (m)^^. The result is in agreement with 
that in FIG. [T] It is worth noting that heavy Majorana 
neutrinos are nearly degenerate in mass, i.e. they form 
a pseudo-Dirac neutrino, so their contributions to {m) 
can be neglected. Interestingly, the next-generation Oz^2/3 
decay experiments are expected to probe (m)^^ with the 
accuracy of 10 — 50 meV, so our model can be tested 
experimentally in the near future. 

c. Search Jar the tri-lepton signals at the LHC. As 
shown in Eq. p2|) . the heavy singlet P couples to the 
gauge sector of the SM, and thus if kinematically acces- 
sible, could be produced at hadron colliders. In the case 
M > Mjj (where denotes the Higgs mass), the heavy 
neutrino can decay in the channels P -I- W~ , P 

D + Z, and P D + H . Now that the heavy neutrinos 
form a pseudo-Dirac particle, we shall focus our atten- 
tion on the Icpton-numbcr-conserving processes induced 
by it. For example, one very interesting and prospec- 
tive channel is the production of three charged leptons 
with missing energy [3l|, i.e. pp — *■ t^t'^t^v{v) -I- jets. 
Another possible one is the pair production of charged 
leptons in different flavors and without missing energy, 
i.e. pp + jets. However, it is difficult to make 

significant observation in this channel at the LHC due to 
the large SM background [3^ . 
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V. SUMMARY 

111 this work, we have considered a novel minimal see- 
saw model, which is quite natural to bring down the 
typical seesaw mechanism to the electroweak scale. Our 
model is minimal among the fermionic seesaw scenarios in 
the sense of particle content, since only two right-handed 
neutrinos are introduced. However, quite different from 
the traditional MSM, the heavy right-handed neutrinos 
in the current consideration is located around the TeV 
scale, and possess sizable Yukawa couplings, which make 
them observable, in particular at the LHC. Furthermore, 
light neutrino masses are protected by a global U{1) sym- 
metry, and thus free from radiative instability. Since the 
model contains only two heavy singlet neutrinos, it is very 
predictive compared to the other seesaw mechanisms. We 
have shown that, in the assumption of a /x r symmetry. 



the tri-bimaximal pattern of lepton flavor mixing can be 
reproduced and only the inverted neutrino mass hierar- 
chy is allowed. The light neutrino masses, the Majorana 
phase and non-unitarity parameters are connected in our 
model, and their non-trivial correlations could be tested 
through a combined analysis of the future long-baseline 
neutrino oscillation experiments, the rare LFV decays, 
the Oi^2/3 decays and the collider signals at the LHC. 



Acknowledgments 

One of the authors (S.Z.) is indebted to G.G. Raffelt 
for warm hospitality at Max-Planck- Institut fiir Physik. 
This work was supported by the Royal Institute of Tech- 
nology (KTH), project no. SII-56510 (H.Z.), and the 
Alexander von Humboldt Foundation (S.Z.). 



[1] P. Minkowski, Phys. Lett. B67, 421 (1977). 

[2] T. Yanagida, in Proc. Workshop on the Baryon Num- 
ber of the Universe and Unified Theories, edited by 
O. Sawada and A. Sugamoto (1979), p. 95. 

[3] M. Gell-Mann, P. Ramond, and R. Slansky, in Supergrav- 
ity, edited by P. van Nieuwenhuizen and D. Freedman 

(1979) , p. 315. 

[4] R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 

44, 912 (1980). 
[5] A. Pilaftsis, Z. Phys. C55, 275 (1992), hep-ph/9901206. 
[6] J. Kersten and A. Y. Smirnov, Phys. Rev. D76, 073005 

(2007), arXiv:0705.3221. 
[7] R. N. Mohapatra and J. W. F. Valle, Phys. Rev. D34, 

1642 (1986). 

[8] M. Malinsky, T. Ohisson, Z.-z. Xing, and H. Zhang, 

Phys. Lett. B679, 242 (2009), arXiv:0905.2889. 
[9] Z.-z. Xing and S. Zhou, Phys. Lett. B679, 249 (2009), 
arXiv:0906.1757. 
[10] Y. Yamanaka, H. Sugawara, and S. Pakvasa, Phys. Rev. 

D25, 1895 (1982). 
[11] H. Zhang, Phys. Lett. B655, 132 (2007), hep- 
ph/0612214. 

[12] P. H. Frampton, S. L. Glashow, and T. Yanagida, Phys. 

Lett. B548, 119 (2002), hep-ph/0208157. 
[13] W.-l. Guo, Z.-z. Xing, and S. Zhou, Int. J. Mod. Phys. 

E16, 1 (2007), hep-ph/0612033. 
[14] S. M. Bilenky and S. T. Petcov, Rev. Mod. Phys. 59, 671 

(1987). 

[15] J. Schechter and J. W. F. Valle, Phys. Rev. D22, 2227 

(1980) . 

[16] G. Lazarides, Q. Shafi, and C. Wetterich, Nucl. Phys. 

B181, 287 (1981). 
[17] R. N. Mohapatra and G. Senjanovic, Phys. Rev. D23, 



165 (1981). 

[18] M. Malinsky, T. Ohisson, and H. Zhang, Phys. Rev. D79, 

011301(R) (2009), arXiv:0811.3346. 
[19] R. Foot, H. Lew, X. G. He, and G. C. Joshi, Z. Phys. 

C44, 441 (1989). 
[20] M. Malinsky, T. Ohisson, and H. Zhang (2009), 

arXiv:0903.1961. 
[21] A. Pilaftsis and T. E. J. Underwood, Phys. Rev. D72, 

113001 (2005), hep-ph/0506107. 
[22] M. Malinsky, J. C. Romao, and J. W. F. Valle, Phys. 

Rev. Lett. 95, 161801 (2005), hep-ph/0506296. 
[23] J. Schechter and J. W. F. Valle, Phys. Rev. D25, 774 

(1982). 

[24] P. F. Harrison, D. H. Perkins, and W. G. Scott, Phys. 

Lett. B530, 167 (2002), hep-ph/0202074. 
[25] Z.-z. Xing, Phys. Lett. B533, 85 (2002), hep- 

ph/0204049. 

[26] G. AltareUi and D. Meloni, Nucl. Phys. B809, 158 

(2009), arXiv:0809.1041. 
[27] T. Schwetz, M. A. Tortola, and J. W. F. Valle, New J. 

Phys. 10, 113011 (2008), arXiv:0808.2016. 
[28] S. Antusch, J. P. Baumann, and E. Fernandez-Martmez, 

Nucl. Phys. B810, 369 (2009), arXiv:0807.1003. 
[29] A. Ilakovac and A. Pilaftsis, Nucl. Phys. B437, 491 

(1995), hep-ph/9403398. 
[30] J. Bernabeu, A. Santamaria, J. Vidal, A. Mendez, and 

J. W. F. Valle, Phys. Lett. B187, 303 (1987). 
[31] F. del Aguila and J. A. Aguilar-Saavedra, Phys. Lett. 

B672, 158 (2009), arXiv:0809.2096. 
[32] F. del Aguila, J. A. Aguilar-Saavedra, and R. Pittau, 

JHEP 10, 047 (2007), hep-ph/0703261. 



